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PROPORTION ESTIMATC PROCESSOR 


1. SCOPE 

This document contains the design of the Proportion Estimate Processor which 
was written to satisfy the software requlreiwnt of Part B of the PI A 
experiment. The purpose of the Proportion Estimate Processor Is to evaluate 
four estimation techniques In order to get an Inproved estimate of the 
proportion of a scene that Is planted In a selected crop. The four techniques 
to be evaluated have been provided by Techniques Development Section and 
are (1) Random Sampling, (2) Proportional Allocation - Relative Count 
Estimate, (3) Proportional Allocation - Bayesian Estimate, and (4) Sequential 
Bayesian Allocation. The user will be given two options for computation of 
the estimated mean square error. These will be referred to as the Cluster 
Calculation option and the Segment Calculation option. 
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2.0 APPLICABLE DOCUMENTS 


2.1 Technical Memorandum: 

The Hulticategory Case of the Sequential Bayesian Pixel Selection and 
Estimation Procedure by M. 0. Pore and T. B. Dennis. 

2.2 LEC- 13945 

Clustering Algorithm Evaluation and the Development of a Replacement for 
Procedure 1 by R. K. Lennington and J. K. Johnson. 

2.3 LEC- 13940 ASA Proceedings 1979 

Bayesian Techniques In Stratified Proportion Estimation by M. D. Pore. 

2.4 LEC-12566 Earth Observation Division Version of 

the Laboratory For Application of Remote Sensing Syst«n (EOD-LARSYS) . 
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3.0 SYSTEMS DESCRIPTION 


3.1 HARDWARE DESCRIPTION 

The software for the Proportion Estimate Processor will be operational on 
the IBM 3031 ccmiputer at Purdue. 

3.2 SYSTEM DESCRIPTION 
3.2.1 OVERVIEW 

The Proportion Estimate Processor will read a control card file (4.2 ), the 
console, and two disk files created by the Pixel Selection and Display 
Processor. These two files (3.2.4. 1) and (3. 2. 4. 2) contain Information 
about the number of clusters, size of each cluster, how the clusters are 
combined, and the labeled dots selected from the clusters. The processor 
will then calculate each of the four proportion estimates required for each 
crop label. The results will be output In a printer report (3. 2. 5. 3) and 
two disk files (3. 2. 5.1) and (3. 2. 5. 2). Appendix A contains a flow of this 
procedure. Fourteen software routines have been developed for this 
processor. Appendix A contains a flow of these routines and others used from 
the EOD LARSYS library routines. 


3.2.2 COMMON BLOCK ESTIM 


The comnon block ESTIM Is used to communicate between the routines In the 
Proportion Estimate Processor. The parameter definitions are as follows: 


NOCLS 

ISUKCL (30) 

NTOT 

NDOTS 


Number of clusters In map file (max * 30). 

Number of pixels In each cluster In the scene. 

Total number of pixels In scene = sum of ISUMCL values. 

Number of dots In the Input label file with labels belonging 
to the selected categories, (number of entries In arrays 
LABEL and ISORT - max 500). 


LABEL(9,500) Information about dots - one entry oer dot. 
(1 ,1) = Analyst Label 
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NLAB(30) 

IS0RT(500) 

NPTS 

P(504,5) 

THRES 

LCNT(5) 

LAB(10,5) 

NOCAT 

ALPHA{5) 

I OPT 


(2.1) > Line nunter 

(3.1) ■ Sample number 

(4.1) > Original cluster number 

(5.1) > Resulting cluster number 

After combining clusters 

(6.1) > Grid dot chosen from 

(7.1) * Indicator for dots used In Random Sample Estimate 

( » 0, not selected; fQ selected) 

(8.1) * Indicator for dots used for Proportion Estimators 

( « 0, not selected; ii<0, selected) 

(9.1) « Indicator for dots used for Sequential Bayesian 

Estimator ( » 0, not selected, ¥0 selected). 

Count of dots In label array for each cluster. 

Sorted pointer table to entries In array LABEL > Sorted by 
ascending cluster number. 

Number of dots user wishes to use In calculating estimates. 

Proportion Estimate Results 

P(1,I) s Random sampling results for label I. 

P(2,I) = Proportional Allocation - Relative Count results for 
label I. 

P(3,I) = Proportional Allocation - Bayesian results for label I. 
P(J,I) = Sequential Bayesian results for label I (J*4 -► KCNT) 
Threshold MSE value to be used In Sequential Bayesian (If zero 
use NPTS to determine estimate). 

Number of labels of Interest In each category. 

Array containing labels of Interest for each category. 

Number of categories for evaluation maximum « 5. 

Weighting factor for each cate^ry to be used In calculating 
mean square errors. 

Options chosen 

> 1 Use a fixed number of dots and use the cluster calculation 
technique. 
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* 2 Use a fixed nutnber of dots and use the segnents calculation 
technique. 

> 3 Use a number of dots determined by the threshold value In 

the Sequential Bayesian and the Cluster Calculation technique. 

> 4 Use a number of dots determined by the threshold value in 

the Sequential Bayesian and the segment calculation techniques. 

FMSE(504) ■ Mean Square Error 

FSME(l) > MSE for Random Sampling 

FSME(2) - MSE for Proportional Allocation Relative Count Estimate. 
FSME{3) ■ MSE for Proportional Allocation Bayesian 
PSME(J), J"4, KCNT « MSE for Sequential Bayesian 

NC0UNT(504) Number of dots used for this calculation. 

NCOUNT (1) ■ NPTS for Random Sampling. 

NCOUNT (2) ■ NPTS for Proportional Allocation Relative Count 
Estimate. 

NCOUNT (3) * NPTS for Proportional Allocation Bayesian. 

NCOUNT (4) = Number of dots used for Sequential Bayesian 
calculation for 2 dots per active cluster. 

NCOUNT (J), J-5. KCNT 

(J»5 NCOUNT (5) » NCOUNT (4) +1 
J«KCNT NCOUNT (J) « NPTS) 

KCNT = Number of entries In NCOUNT, 

FMSE, and P arrays. 

A(5) « Input values for constant A for each category - used In Bayesian 

calculation - required Input for more than 2 categories. 

3.2.3 PROGRAM DOCUMENTATION 

3.2. 3.1 PRPREST 

Purpose 

This Is the driver routine for the Proportion Estimate Processor. 
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Linkages 

PRPEST calls SETEST. RANDOM. PROPOR. and BAYES. 


Interfaces 

Interface Is accomplished through coimon block ESTIM and subroutine calling 
arguments . 

Inputs 

Outputs 


Storage Requirements 


Description 

The PRPEST program coordinates the logical steps In calculating the 
proportion estimate. It calls SETEST to read In the control cards and 
data files. It calls RANDOM to calculate the Random Sampling Proportion Estimate. 
It calls BAYES to calculate the Sequential Bayesian Allocation Proportion 
Estimate. It calls PROPOR to calculate both the Relative Count Estimate 
and the Bayesian Estimate for Proportional Allocation. It writes both 
bhe output Dot File and Results File and writes an EXEC file for passing 
fUetype names back to the controlling EXEC file. 

Flow Chart 
Reference listing. 

Listing 

See Appendix B.l . 
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3.2. V2 SETEST 


Purpose 

The purpose of the SETEST subroutine Is to read the control cards and the 
tMO disk files (3. 2. 4.1) and (3. 2. 4. 2) and start to build the conmon biiock 
ESTIM. 


linkages 

SETEST Is called by PUREST and It calls the EODLARSYS subroutines NXTCHR, 
NUWER. and FLTNUM 

Interfaces 

SETEST Interfaces with other rout nes by use of calling arguments and the 
common block ESTIM. It also uses the common block QIOBM as a system 
standard Interface. 


Inputs 


Calling sequence: CALL SETEST (lERR, ISEG, ITYPE, IXCNT, IXLAB, ICLCNT, ICLS) 


PARA«TER 

DIMENSION 

IN/OUT 

DEFINITION 

lERR 

1 

OUT 

Error Indicator for errors In set up 
processing. 

■ 0 no error 

■ 1 error encountered 

ISEG 

1 

OUT 

Se^nent number for which Estimate Is 
performed. 

ITYPE 

1 

OUT 

Type of dots used NXXY. 
N ■ R * reformatted 
N ■ G ■ ground truth 
N ■ I ■ Integrated 
XX ■ analyst Initials 
Y ■ version number 

IXCNT 

1 

OUT 

Count of labels to be Ignored by 
processor range 1 to 3. 

IXLAB 

3 

OUT 

Labels to be Ignored (the label X will 
always be Ignored J 
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PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

ICLCNT 

1 

OUT 

Count of number of original Input 
clusters that have been combined 




Into other clusters. 

ICLS 

30 

OUT 

Cluster numbers that have been 
combined Into other clusters. 


The SETEST reads the control cards documented In section 4.2. 

The SETEST read ^e two disk files documented In section 3.2.4. 1 and 3. 2. 4. 2. 
SETEST reads the console for the type of dots used. 

Outputs 

Report of Input cards summary and error messages. 

Storage Requirements 


Description 

The SETEST routine reads the control cards and the tw Input data files. 

It sets up the rer^<id buffer and reads the keyword and data fr(»n the cards 
checking for errors and printing error messages as applicable. The 
subroutine then reads the cluster Information data file filling In NOCLS 
and ISUMCL, & NTOT In ESTIM. It then checks for combined clusters and 
reworks ISU^L If necessary. If tM> clusters are combined, the remaining 
cluster’s ISUMCL entry will be the sum of the entries for the two combined 
cluster, and the eliminated cluster ISUMCL entry will be zero. Next the 
labelled dot data file will be read Into the array label. If a label Is 
blank. X. or matches a label Input on the IGNORE card, or does not match an 
Input label, then that dot will be Ignored. The number of dots accepted for 
the label array will become NOOTS. Next the ISORT array will be bulU to 
contain printers to the dots In the label array by ascending cluster number 
order. The arr^y IRES will be checked for combined clusters so that the 
pointers can be grouped to Include cluster numbers for combined clusters. 
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While building the ISORT array the NLAB array of count of dots for each 
cluster remaining after combining will be built. The remaining clusters 
will be referred to as active clusters. 


Flow Chart 
Reforence listing. 

Listing 

See Appendix B.2. 
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3. 2. 3. 3 BAYES 


PurDOse 

The purpose of the BAyos subroutine Is to calculate the Sequential Bayesian 
Allocation Proportion Estimate. 


Linkages 

BAYES Is called by PRPEST and calls the random number generator subroutine 
RANDU and the subroutines FMSES, FMSEC, DMSES, and DMSEC. 

Interfcces 

BAYES Interfaces with other routines by use of calling arguments and the 
common block ESTIN. 


Inputs 


Calling sequence: CALL BAYES (IPICK, 

IX, NI, IXI) 

PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

IPICK 

500 

IN/OUT 

Work array to be used for picking 
random pointers Into ISORT array 
which then point to entries In the 
LABEL array. 

IX 

1 

IN/OUT 

Seed value for subroutine RANDU 
changed after each call to RANDU. 

NI 

30 

IN/OUT 

Work array to contain counts of 
dots to be chosen from each active 
cluster. 

IXt 

30,5 

IN/ OUT 

Work array to contain count of dots 
chosen in each cluster which have 
a label of Interest for each 
category. 


Output 

Printer report of mean square error and proportion estimate values as points 
are picked. The report will also specify which cluster each additional dot 
Is picked from. 
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s torage Requirements 


Description , 

The subroutine BAYES first randbmiy picks two dots in each active cluster. 

If there are only two categories of interest, then the subroutine will 
recalculate the values for the constants A used in the calculations by 
calculating a PINIT with the formula PINIT(I) - - (ISUMCL(J)/NOT)* 
{IXI(J,I)+A(I)+1)/(NI{J)+ASUM + NOCAT) >1 

where ASUM ■ A(I) 

If PINIT(l) greater than .5 

If PINIT(l) less than ,5 

JpiNITd) 1-1 , A(2) • 0 

. “[i - MNTrri)] 

When more than two categories are involved, the user must have input the 
A(I) to use. Now the subroutine enters a processing loop to do the following: 

1. Calculate the mean square error by calling FKSES or FMSEC depending on 
the calculation option chose. 

2. Calculate the proportion estimate as: 

ISUMCUJ) * (IXI(J.I) + A(I) +1 

liTcr Nnjy+Aiuff+ n6 

for each category I. 

3. Check that all necessary dots have been picked either by fixed number 
of dots or by MSE threshold reached. If all dots have been picked, 
return to PRPREST. 

4. If more dots are needed calculate the AMSE value for each cluster by 
calling either DMSES or DMSEC depending on the calculation option chosen. 

5. Determine which cluster to pick the next dot from by deciding which 
cluster has the largest AMSE value. 

6. Pick a dot from the required cluster recalculate IXI and NI and return to 
point 1 in the loop. 



LEKSCO- 15353 


3-9 


Flow Chart 


fteference listing. 
listing 

See Appendix B.3. 
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3.2. 3.4 RANDOM 


Purpose 

The purpose of the RANDOM subroutine Is to compute the Random Sample 
Proportion Estimate. 

linkages 

RANDOM Is called by PRPEST and It calls the random number generator 
subroutine RANDU. 

Interfaces 

RANDOM Interfaces with other routines by use of calling arguments and the 
coimion block ESTIM. 

Inputs 

Calling sequence; CALL RANDOM (IPICK.IX) 


PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

IPICK 

500 

IN/OUT 

Work array to be used for picking 
random pointers Into LABEL array. 

IX 

1 

IN/ OUT 

Seed value for subroutine RANDU 
changed for each call to RANDU. 

Outputs 





Report of results of calculation for Random Sampling Proportion Estimate. 
Storage Requirement 


Description 

The subroutine RANDOM calls the subroutine RANDU to generate random numbers, 
RANDOM scales them to a dot number entry In If^Zl and verifies that each dot 
Is only picked one time. The RANDOM counts the number of chosen dots (NPTS) 
LEMSCO- 15353 


with a label I of Interest (XLAB) and calculates the proportion estimate P(1» 
where P(1,I) ■ XLAB/NPTS for category I. Then the mean square error Is 
calculated as FMSE(l) - ALPHA(I)*(P(1 ,1)*(1-P(1 .D)/(NPTS-1) 


I-l 


I) 


Flow Chart 
Reference listing. 


Listing 

See Appendix B. 4. 
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Purpose 


The purpose of the PRW*0R subroutine Is to calculate the Proportional 
Allocation - Relative Count Estimate and the Proportional Allocation Bayesian 
Estimate. 

Linkages 

PROPOR is called by PRPEST and calls the subroutines RANDU, FMSEC, and 
FMSES. 

Interfaces 

PROPOR interfaces with other routines by use of calling arguments and the 
common block ESTIM. 

Inputs 


Calling Sequence; CALL PROPOR(IPICK.IX.NI.IXI) 


PARAMETER 

DIMENSION 

IN/ OUT 

DEFINITION 

IPICK 

500 

IN/OUT 

Work array to be used for picking 
random pointers into ISORT array 
which then points to entries in 
the LABEL array. 

IX 

1 

IN/OUT 

Such value for subroutine RANDU - 
changed for each call to RANDU. 

NI 

30 

IN/OUT 

Work array to contain count of 
dots to be chosen from each 
active cluster. 

IXI 

30,5 

IN/ OUT 

Work array to contain count of 
dots chosen in each cluster which 
have a label of interest for each 
category. 


Outputs 

Report of results of calculations for Proportional Allocation Relative Count 
Estimate and for Proportional Allocation Bayesian Estimate. 
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storage Requirements 


Description 

PROPOR first calculates the number of dots to pick from each active cluster. 
It uses the formula NI(J) = NPTS * ISUMCL (J)/NT0T where these variables 
are defined in section 3.2.2. Next the sum of NI(J) is checked to be sure 
it equals NPTS. If not, when the NI values are adjusted up or down by one 
starting with the largest NI value until NI(J) * NPTS. Next NI dots 

are chosen for each cluster J; they are chosen randomly and checked for no 
duplication. The labels of interest are counted in array IXI (J,I). Next 
the Proportional Allocation Relative Count Estimate and its n?ean square error 
are claculated by the following formulas: 

and 

FMSE(2) = ALPHA(l)* / p(2,I) * (1-P(2,I)) ) 

Z- V NPTS-1 / 

1=1 

Then for the two category case calculate values for A as 


for Pf2,l) less than .5 for P(2,l) greater than .5 

A(,) . = » *(') = »• 


Use the user input A values for more than two categories. Using the A 
values, calculate the Proportional Allocation Bayesian Estimate 


P(3,I) 



NOCAT 


flSUMCL(J) ] 

*/ 

\^T0T / 

\ 


IXI(J,I) A(I 
NI ( J ) + ASUM+N 



where ASUM = A(I). Finally calculate the Bayesian mean 


by calling either FMSES or FMSEC according to the calculation 


square error 
option chosen. 
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Flow Chart 


Reference listing. 
Listing 

See Appendix B.5. 
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3.2. 3.6 RANDU 


Purpose 

The subroutine RANDU will generate random numbers for use by the Proportion 
Estimator routines. 

Linkages 

RANDU Is called by RANDOM. PROPOR. and BAYES. 

Interface 

RANDU Interfaces with other routines by use of calling arguments. 

Inputs 

Calling Sequence: CALL RANDU(IX.IY.YFL) 


PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

IX 

1 

IN 

Seed value for generating random 
number. 

lY 

1 

OUT 

Integer value to be used for IX 
on next call to RANDU. 

YFL 

1 

OUT 

Random number between 0 and 1 
generated by RANDU. 


Storage Requirements 


Description 

RANDU takes the seed value IX and multiplies It by 65539. Then It converts 
to a positive Integer If necessary. Finally It scales the Integer to a 
real number between 0 and 1. A suggested start value for the IX seed value 
Is 187521429 but another nine digit Integer can be used. 
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Flow Chart 
Reference listing. 


Listing 

See Appendix B.6. 
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3. 2. 3. 7 FMSES 


Purpose 

The subroutine FMSES will calculate the Keans Square Error for the segment 
calculation option for the Bayesian estimates. 

Linkages 

FMSES Is called by PROPOR and BAYES and It calls the function subroutines 
THETA, BIAS, and VAR. 

Interfaces 

FMSES Interfaces with other routines through calling arguments. 

Inputs 

Calling Sequerce: CALL FMSES (FMSE, ALPHA. IXI. Nl, A, ASUM. NOCAT. NOCLS. 

ISUMCL, RTOT) 


PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

FMSE 

1 

OUT 

Mean Square Error calculation 
result. 

ALPHA 

5 

IN 

Weighting factor for each category 
of Interest. 

IXI 

30,5 

IN 

Count of dots chosen for each 
category from each cluster. 

NI 

30 

IN 

Count of dots chosen from each 
cluster. 

A 

5 

IN 

Constant A for use In MSE equation 
one value per category. 

ASUM 

1 

IN 

Sum of A values. 

HCCAT 

1 

IN 

Number of categories of interest. 

NOCLS 

1 

IN 

Number of clusters In segment. 

ISUMCL 

30 

IN 

Count of pixels In each active 
cluster In segment. 
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PARAHETER 


DIMENSION IN/OUT 


DEFINITION 


RTOT 


IN Total number of pixels In segment 

• sum of ISUMCL values. 


Output 

Retults are returned In calling argument. 
Storage Reoulrwients 


Description 

This subroutine Mill calculate the mean square error according to the 
equation 

NOm 

FMSE « L alpha (I) * CMSEd ) 

I»1 


where: 

CMSE(I) - V(I)+B(I)**2 


V(I) 

and 


■ r [( 


mis t 

B(I) * E ^ ISMUCL(J) Y BIAS(J) 


ISUMCL(J) **2 

■”RTor^ 


)]■ 


VAR(J) 


It uses the function subroutines THETA, VAR, ind BIAS to get the value THETA 
and the values for VAR(J) and BIAS(J) used In the equations above. 


Flow Chart 
Reference listing. 


Listing 

See Appendix B.7. 
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3. 2. 3. 8 FKSEC 
Purpose 

The subroutine FKSEC will calculate the Kean Square Error for the Cluster 
Calculation option for the Bayesian Estimators. 

Linkages 

FI^EC Is called by PROPOR and BAYES and It calls the function subrriutines 
THETA and RKSE. 

Interfaces 

FKSEC Interfaces with the other routines through calllmg arguments. 

Inputs 

Calling Sequence: CALL FKSEC (FKSE. ALPHA, IXI, NI, A, ASUK. NOCAT. NOCLS, 

ISUKCL, RTOT) 


PARAMETER 

DIKENSION 

IN/OUT 

DEFINITIONS 

FMSt 

1 

OUT 

Mean Square Error calculation 
result. 

ALPHA 

5 

IN 

Weighting factor for each category 
of Interest. 

IXI 

30,5 

IN 

Count of dots chosen for each 
category from each cluster. 

NI 

30 

IN 

Count of dots chosen for each 
cluster. 

A 

5 

IN 

Constant A for use In MSE equation 
one value per category. 

ASUK 

1 

IN 

Sum of the A values. 

NOCAT 

1 

IN 

Number of categories of interest. 

NOCLS 

1 

IN 

Ntanber of clusters In se^nt. 

ISUKCL 

30 

IN 

Count of pixels In each cluster In 
se^nt. 
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PARAMETER DIMENSION IN/OUT 


DEFINITION 


RTOT 1 IN Total number of pixels in segment 

« sum of ISUMCL values. 


Outputs 

Results are retained In calling argument. 
Storage Requirements 


Description 


This subroutine will calculate the mean square error according to the 
equation FMSE ALPHA(I) * CMSE(i) 


E 

where 
CMSE(I) = 


NncLS 

E 


niSUMCL(J)l**2\*RMSE(J) 
H RTOT I 


The subroutine was the function subroutine RMSE to calculate the value for 
RMSE(J) in the equation above. It use the subroutine THETA to calculate 
the value for THETA to Input to RMSE. 


Flow Chart 
Reference listing. 


Listing 

See Appendix B.8. 
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3.2.3.9 OMSES 
Purpose 

The purpose of OMSES Is to calculate the value of DMSE for each cluster as 
requested. When the segment calculation option Is chosen for the Sequential 
Bayesian Allocation. 

Linkages 

OMSES Is called by BAYES and calls the function subroutines THETA, VAR, and 
BIAS. 

Interfaces 

OMSES Interfaces with other routines through calling arguments. 

Inputs 

Calling Sequence: CALL OMSES (IIRST, ILST, ISUMCL, IXI, NI, ASUM, NOCAT, A, 




ALPHA, 

RTOT 

, OMSE) 

PARAMETER 

DIMENSION 

IN/OUT 


DEFINITION 

IFRST 

1 

IN 


First cluster to calculate AMSE 
value for. 

ILST 

1 

IN 


Last cluster to calculate AMSE 
value for. 

ISUMCL 

30 

IN 


Count of pixels In each clusters 
In segment. 

IXI 

30,5 

IN 


Count of dots chose for each 
category from each cluster. 

NI 

30 

IN 


Count of dots chosen for each 
cluster. 

ASUM 

1 

IN 


Sum of A values. 

NOCAT 

1 

IN 


Number of categories of interest. 

A 
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IN 
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Constant A for use In AMSE 
calculations 1 per category. 


PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

ALPHA 

5 

IN 

Weighting factor for each category 
of Interest. 

RTOT 

1 

IN 

Total number of pixels In segment 
■ sum of ISUMCL values. 

DHSE 

30 

OUT 

DMSE value calculated for each cluster 
requested. 

Outputs 





The results are returned by the calling argument. 


Storage Requirements 


Description 

This subroutine will calculate a AMSE value for each cluster from IFRST to 

ILST according to the following equation: 

NOCAT 

DMSE(J) = 2 ALPHA(I) * CMSE(J.I) 

' I-l 

where 

CMSE(J.I) - RMSEl - (d-Tl) ♦ RMSE2) - T1*RMSE3 
RMSE 1 « MSE (THETA(NI(J), IXI(J,I))) 

RMSE2 » MSE (THETA(NI(J)+1 , IXI(J,I))) 

RMSE3 » MSE (THETA(NI{J)+1 , IXI(J,I)+I)) 

T1 = THETA (NI(J), 1X1(0, 1)) 

MSE « V + B**2 

V * 

B. *BIAS(J) 

OMSES uses the function subroutines THETA, VAR, and BIAS(J), and THETA(N,X) 
used above. 
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Flow Chart 
fteference listing. 

Listing 

See Appendix B.9. 
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3.2.3.10 DMSEC 
Purpose 

The purpose of DMSEC Is to calculate the value of AMSE for each cluster as 
requested when the cluster calculation option Is chosen for the Sequential 
Bayesian Allocation. 

Linkages 

DMSEC Is called by BAYES and calls the function subroutines THETA and RMSE. 
Interface 

DMSEC Interfaces with other routines through calling arguments. 

Inputs 

Calling Sequence: CALL DMSEC(IFRST. ILIST, ISUMCL, IXI, NI, ASUM, NOCAT. A, 

ALPHA, RTOT, DMSE) 


PARAMETER 

DIMENSION 

IN/ OUT 

DEFINITION 

IFRST 

1 

IN 

F’rst cluster to calculate AMSE 
value for. 

ILST 

1 

IN 

Last cluster to calculate AMSE 
value for. 

ISUMCL 

30 

IN 

Count of pixels In each cluster 
for segment. 

IXI 

30,5 

IN 

Count of dots chosen for each 
category from each cluster. 

NI 

30 

IN 

Count of dots chosen for each 
cluster. 

ASUM 

1 

IN 

Sum of A values. 

NOCAT 

1 

IN 

Number of categories of Interest. 

A 

5 

IN 

Constant A for use In AMSE 
calculations - 1 per cluster. 

ALPHA 

5 

IN 

Welghtine factor for each category 
of Interest. 
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PARAMETER 

DIMENSTION 

I M/OUT 

DEFINITION 

RTOT 

1 

IN 

Total number of pixels in segment 
- sum of ISUMCL values. 

DMSE 

30 

OUT 

AMSE values calculated for each 
cluster requested. 

Outputs 





The results are returned by the calling argument. 
Storage Requirements 


Description 

This subroutine will calculate aAKSE value for each cluster from IFRST to 
ILST according to the following equation: 

DMSt(J) = ALPHA(I) * CMSE(J.I) 

u 1=1 

where p -| 

CMSE(J.I) = (^ SUMCL(J) **2)l * (RMSF.l - (1-Tl)* RMSE2 - T1*RKSE3) 

RMSE1 = RMSE (THETA(NI(J),IXI(J,I))) 

RMSE2 = RMSE (THETA(NI(J)+1 , IXI(J,I))) 

RMSE3 = RMSE (THETA(NI{J)+1 ,IXI(J,I)+1 )) 

TI = THETA(NI(J),IXI(J,D) 

DMSEC use the function subroutines THETA and RMSE to calculate the values 
above. 

Flow Chart 

Reference listing. 

Listing 

See Appendix B.10. 
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3.2.3.11 THETA 
Purpose 

The purpose of the function subroutine THETA Is to calculate the value of 
THETA to be used for computation of mean square errors and AMSE In the 
Bayesian Estimation Segment Calculation option. 

Linkages 

THETA Is called by subroutines FMSES, FMSEC, DMSES, and DMSEC. 

Interface 

THETA Interfaces with other routines through calling arguments. 


Inputs 

Calling Sequence: THETA(IXI,NI.A,ASUM,NOCAT) 


PARAMETER 

DIMENSION 

IN/OUT 

DEFINITION 

IXI 

1 

IN 

Count of number of dots chosen 
from NI dots with label of Interest. 

NI 

1 

IN 

Count of total nun^er of dots 
chosen for this cluster. 

A 

1 

IN 

Input constant A for this cateogry. 

ASUM 

1 

IN 

Input constant which Is sum of A 
for all . 

NOCAT 

1 

IN 

Number of categories of Interest. 


Outputs 

The value of THETA is returned to the calling routine. 




Storage Requirements 




1 
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Description 


The function subroutine THETA will calculate the value for THETA according 
to the following equation: 


THETA » 


NI+ASUWfNOCAT 


Flow Chart 
Reference listing. 


Listing 

See Appendix 6.11 . 
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3.2.3.12 BIAS 


Purpose 

The purpose of the function subroutine BIAS Is to calculate the value of 
the bias for use In calculations of mean square error and OMSE for the 
Bayesian Estimators Segment Calculation option. 

Linkages 

BIAS Is called by FMSES and DMSES. 

Interfaces 

BIAS Interfaces with other routines through calling arguments. 

Inputs 

Calling Sequence: BIAS(THETA.NI .A.ASUM, NOCAT) 


PARAMETER 

DIMENSION 

IN/OUT 

DESCRIPTION 

THETA 

1 

IN 

Value for THETA calculated by 
routine THETA. 

NI 

1 

IN 

Count of total number of dots 
chosen for this cluster. 

A 

1 

IN 

Input constant for this category. 

ASUM 

1 

IN 

Sum of A values used for all 
categories. 

NOCAT 

1 

IN 

Number of categories of Interest. 


Outputs 

The value of BIAS Is returned to the calling subroutine. 
Storage Requirements 
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Description 

The function subroutine BIAS Mill calculate the value of the bias with 
the following equation: 

. A- M - THETA * (ASUiHNOCAT) 

NI+ASUMfNOCAT 


Flow Chart 
Reference listing. 


Listing 

See Appendix B.12. 
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3.2.3.13 VAR 
Purpose 


The purpose of the function subroutine VAR Is to calculate a value for 
variance for use In calculations of mean square error and AMSE for the 
Bayesian Estimators Segment Calculation option. 

Linkages 

VAR Is called by FMSES and DMSES. 

Interfaces 

VAR Interfaces with other routines through calling argument. 

Inputs 

Calling Sequence: VAR(THETA.NI.A.ASUM.NOCAT) 


PARAMETER 

DIMENSION 

IN/OUT 

DESCRIPTION 

THETA 

1 

IN 

Value for THETA calculated by 
routine THETA. 

NI 

1 

IN 

Count of total number of dots 
chosen for this cluster. 

A 

1 

IN 

Input constant for this category. 

ASUM 

1 

IN 

Sum of A values used for all 
categories. 

NOCAT 

1 

IN 

Number of categories of Interest. 


Output 

The value of VAR h returned to the calling subroutine. 
Storage Requirements 
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Description 


The function subroutine VAR will calculate the value of the variance with 
the following equation: 



Flow Chart 
Reference listing. 

Listing 

See Appendix B.13. 
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3.2,3.14 RMSE 
Purpose 

The purpose of the function subroutine RMSE Is to calculate the Mean Square 
Error for the Sequential Bayesian Cluster Option. 

Linkages 

RMSE Is called by FMSEC and DMSEC. 

Interfaces 

R^E Interfaces with other routines through calling arguments. 

Inputs 

Calling Sequence: RMSE(THETA.IXI .NI,A.ASUH,NOCAT) 


PARAMETER 

DIMENSION 

IN/OUT 

DESCRIPTION 

THETA 

1 

IN 

Value for THETA calculated by 
routine THETA. 

IXI 

1 

IN 

Count of set of NI dots In category 
of Interest. 

NI 

1 

IN 

Count of total number of dots chosen 
for this cluster. 

A 

1 

IN 

Input constant for this category. 

ASUM 

1 

IN 

Sum of A values used for all 
categories. 

NOCAT 

1 

IN 

Number of categories of Interest. 



The value for RMSE Is returned to the calling routine. 


Storage Requirements 
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Description 

The function subroutine RMSE will calculate the value of the mean square 
error according to the following equation: 

NI*THETA*n -THETA) ♦ (A+1 -THETA* (ASUM+NOCAT))^ 

(Nl+A5UMiNGCAT)i 

Flow Chart 
Reference listing. 

Listing 

See Appendix B.14. 
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3.2.4 INPUT FILE FORMATS 


3. 2.4.1 Label File 

The label file Is a card Image file containing a header card, one data card 
per dot labeled, an *END card, and some grid Information cards. Only the 
HEADER, DATA, and *END cards are of Interest to the processor and their 
formats follow. The filename of this file will be SE6M PIANXXY where 
SECW = Segment Number N*R,6 or I XX * Analyst Initials and Y ■ version 
number. 


HEADER CARD 


COLUMN 

FORMAT 

CONTENTS 

1-4 

14 

Number of Data 

5-30 

26A1 

'Pixels select( 

31-35 

14 

Segment Njmber 

DATA CARD 

COLUMN 

FORMAT 

CONTENTS 

1-4 

1A4 

Analyst Label 

5-8 

H4 

Line Number 

9-12 

114 

Sample Number 

13-16 

114 

Cluster Number 

17-20 

114 

Grid Number 

♦END CARD 

COLUMN 

FORMAT 

CONTENTS 

1-4 

1A4 

♦END 
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3. 2. 4. 2 Cluster Information File 


The cluster information file Is a card Image file. The filename of this 
file Is SE6M PR0C2 where SEGM « Segment number. 


CARPI 

COLUMN 

FORMAT 

CONTENTS 

1-4 

15 

Number of cluster In file 

CARD 2 

COLUMN 

FORMAT 

CONTENTS 

1-60 

15A5 

Number of pixels In clusters 1-15. 

CARD 3 

COLUMN 

FORMAT 

CONTENTS 

1-60 

15A5 

Number of pixels In cl users 16-30 

CARD 4 

COLUMN 

FORMAT 

CONTENTS 

1-60 

15A6 

Resulting cluster number after combining for 
clusters 1-15. 

CARD 5 

COLUMN 

FORMAT 

CONTENTS 

1-60 

15A5 

Resulting cluster number after combining for 
cluster 16-30. 
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3.2.5 OUTPUT FILE FORMATS 
3. 2. 5.1 Output Dot File 

The output Dot File is a card image file containing a header card and one 
card per labeled dot used by the Proportion Estimate Processor. The 
filename for this file is SEGM DAANXXY where AA - user options FC, FS, 
MC. or MS and NXXY » the NXXY from the Label File (3. 2. 4.1) 


HEADER CARD 
COLUMN FORMAT 

1-4 1A4 

5 IX 

6 lAl 

7-8 1A2 


9-12 1A4 


CONTENTS 
Ses^nt Number 
Blank 
‘D* 

Processing options used 
'FC for fixed cluster option 
'FS' for fixed segment option 

'MC for MSE cluster option 
'MS' for MSE segment option 

ITYPE - NXXY 

N » type of dots used R, G, or I 
XX » analyst initials 
Y » version number 


13 

IX 

blank 

14-16 

13 

Number of dots used * NPTS 

17 

lY 

Blank 

18-19 

12 

Number of categories 
« NOCAT 

20 

IX 

Blank 

21,28, 

35,42,49 

lAl 

Category label name 
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coLum 

FORMAT 

CONTENTS 

22.29,36. 

lAl 

1 1 

43.50 



23-26, 

30-33, 

37-40, 

44-47, 

51-54 

4A1 

Category labels grouped into category of Interest label 

27.34, 

41,48. 

55 

IX 

Blank 

56-58 

3A1 

Labels to Ignore 

59 

IX 

Blank 

60-61 

12 

Number of clu combined into other cluster 

62 

IX 

Blank 

63,-64, 

66-67, 

69-70, 

72-73, 

75-76 

12 

Cluster numbers combined into other cluster 

65.68,71. 

74,77 

IX 

Blank 
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MTA CARD 


COLUMN 

FORMAT 

CONTENTS 

1-4 

1A4 

Analyst Label 

5-8 

114 

Line Number 

9-12 

114 

Sample Number 

13-16 

114 

Original cluster number 

17-20 

114 

Resulting Cluster number after combining 

21-24 

114 

Grid number 

25-28 

114 

Indicator for dot used In Random Sample Estimate 
= 0 not selected; #0, selected 

29-32 

114 

Indicator for dot used in Proportion Estimates 
= 0 not selected; #0, selected 

33-36 

114 

Indicator for dot used In Bayesian Estimate 
= 0 not selected; #0, selected 



3. 2. 5. 2 Results File 


The Results fl'le Is a card Image file containing one header card and 
several data cards. The first data card will contain the Random Sampling 
Technique results. The second data card will contain the Proportional 
Allocation - Relative Count results. The third data card will contain 
the Proportional Allocation > Bayesian results. The remaining data cards 
will contain the Sequential Bayesian results. The fourth data card will 
contain the results obtained after picking two dots per cluster, the fifth 
will contain the results after adding one dot, etc, until either the 
Threshold MSE value Is reached or the fixed number of dots have been 
chosen. The name of this file will be SEGM RAANXXY where AANXXY matches 
the ones for the output dot file (3. 2. 5.1). The card formats follow. 


Header: 

The header for this file Is the same as the header for the output dot file 
except that column 6 will contain an R for results Instead of the D for 
dots. 


Data; 


COLUW 

FORMAT 

CONTENTS 

1-3 

13 

Number of dots used to compute this estimate. 

4-11 

F8.4 

MSE value 

12-19 

F8.4 

Estimate for category 1 

20-27 

F8.4 

Estimate for category 2 

28-35 

F8.4 

Estimate for category 3 

36-43 

F8.4 

Estimate for category 4 

44-51 

F8.4 

Estimate for category 5 
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3. 2. 5. 3 Output Printer Report 

The output printer report will be in tno parts. Part 1 will be a summary 
of the input control cards and error messages or warning messages as 
required. Part 2 will contain the results of the run in the following order: 

(1) Sequential Allocation Bayesian Estimate) (2) Random Sampling Estimate, 
(3) Proportional Allocation - Relative Count Estimate, and (4) Proportional 
Allocation Bayesian Estimate. See Appendix A for example report. 
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4.0 Users Guide 


4.1 EXEC FILE 

The user will access the Proportion Estimate Processor by use of the PRP 
exec (Appendix B.15). 

To run the processor for a given segment and label file the user will enter 
PRP 1234 PIANXXY 
Where 1234 = segment nurrber 

PIANXXY = Input Label File flletype 
N « Dot set used R,G, or I 
XX « Analyst Initials 
Y « Version nunber 

The following files will be expected on the users A disk. 

1234 PIANXXY = Input Label File (section 3.2.4. 1) 

1234 PR0C2 * Cluster Information File (section 3. 2. 4. 2) 

PRP CC = Control Card File (section 4.2) 

Upon completion of the run the user should find a report output to the 
HOUSTON printer and two output disk files on the A disk: 

1234 DAANXXY = Output Dot File (section 3. 2. 5.1) 

1234 RAANXXY = Output Results File (section 3. 2. 5. 2) 
where AA * User Options Chosen: 

FC.FS.MC. or MS 


4.2 CONTROL CARDS 

The control cards will contain a keyword in columns 1-10 and data starting 
after column 11. 

One of the following two cards is required: 
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Keyword 

NPTS 


THRES 


Default 

none 


none 


Integer 


Real 


Data 

Number of dtos to use for calculation 
of estimates - If this card is chosen, 
the user is chosing the fixed dot 
procedure option and the Sequential 
Bayesian will add dots until this 
number is reached and then stop. 

Threshold MSE value for the Sequential 
Bayesian - if this card is chosen, 
the user has chosen the MSE option and 
the Sequential Bayesian will add dots 
until the MSE value drops below this 
threshold value - when that point is 
reached the number of dots (NPTS) to 
use for the remaining techniques will 
be fixed at the number of dots used for 
the Sequential Bayesian. 


The following Cards are Required 


Keyword 

ALPHA 


Default 


none 


Jm. 

Real 


LABEL 


none 


Character 


Data 

Weighting value for categories to be 
used in calculating MSE values > one 
value per category. 

Labels of categories to be used for 
estimate calculations. The data will 
be in the form 
C*A,BB;N;S 

This says that these are three categories 
of interest: C,N, and S. The category 
C is made of dots labelled C,A, and BB. 
Two character labels are allowed and 
each category (a maximum of u/ ^.an be 
made of up to 10 labels. 
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Keyword Defaul t 


Type Data 

*END none Character No data on this card the keyword signals 

the end of the control cards for this 
run. 


The following card Is 

Keyword Defaul t 

A for 2 

cateog ry 

by one value per category see sections 
3. 2. 3.1 through 3.2.3.14 
for the A usage In equations 


required for more than two categories to evaluate. 
Type Data 

Real The constants A for use In Bayesian 

Technique computations ■ there should 


The following cards are optional 


Keyword 

Default 

Jm. 

Data 

OPTION 

C 

Character 

Calculation technique option. 
C * Cluster Option 
S = Segment Option 

IGNOPI 

none 

Character 

Labels to Ignore In the label file - 
the label X will always be Ignored - 
the user may choose two auditlonal 
labels to Ignore. 

COMMENT 

none 

Character 

Any conment to be used In report heading. 

HEDl 

none 

Character 

Header line 1 for report. 

HE 02 

none 

Character 

Header line 2 for report. 

DATE 

none 

Character 

Date for report. 
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